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Abstract 

Paradoxes in the Boltzmann kinetic theory are presented. Firstly, it 
is pointed out that the usual notion concerning the perfect continuity of 
distribution function is not generally valid; in many important situations 
using certain types of discontinuous distribution functions is an absolute 
must. Secondly, it is revealed that there is no time reversibility in terms 
of beam-to-beam collisions and, in connection with this, there are in- 
trinsic difficulties in formulating the net change of molecular density due 
to collisions, either in the three-dimensional velocity space or in the six- 
dimensional phase space. With help of simple examples, the paradoxes 
manifest themselves clearly. 



1 Introduction 

When the Boltzmann equation, thought of as the first equation of kinetic theory, 
came out, many philosophical doubts arose[l], of which most were related to the 
fact that while Newton's equations themselves were time-reversible the Boltz- 
mann equation, as a consequence of Newton's equations, was time-irreversible. 
As time went by, in particular after the formally rigorous BBGKY theory was 
formulated in the middle of the last centrary[2— 8], concerns of the time-reversal 
paradox gradually faded out, and it was believed that the ultimate understand- 
ing of the related issues had been completely established (at least in the regime 
of classical physics). 

Nevertheless, later developments of mathematics and physics have suggested 
that things should be more complicated than they appear to be. As one of the 
developments, the studies of fractals[9, 10] illustrate that complex structures on 
very small scales resembling those on large scales exist quite generally. Impact 
of the fractal studies on our customary thought can be various and enormous: 
some of us may be led to think about the possibility that nature gets certain 
mechanisms which allow macroscopic structures to shape and form microscopic 
structures, and others may be tempted to construct new formalisms describing 
dynamics of complex structures over infinitely small ranges. 
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As another of the developments, it has been revealed that the standard 
treatment of collisions involves essential problems and the collective behavior of 
molecular collisions should be reinvestigated [11]. 

In this paper, we wish to present a rather comprehensive description of para- 
doxical aspects of the Boltzmann equation. Firstly, it is shown that many mech- 
anisms in a real gas are capable of producing abnormal molecules not associated 
with continuous distribution functions; as long as macroscopic noncquilibrium 
exists somewhere, sharp, even discontinuous, microscopic structures will result 
elsewhere and certain types of discontinuous distribution functions need to be 
employed. Secondly, it is illustrated that formulating the number of molecules 
entering and leaving a volume element is a very tricky task to accomplish. Much 
contrary to the customary thought, there is no time reversibility that we can 
utilize. Even worse, if we want to know how many molecules enter and leave 
an infinitesimal phase volume element, we shall meet with conflicting differ- 
ential processes. All these paradoxes point to one thing: it is not possible to 
have a formalism that is mainly based on differential-type analyses and can still 
describe kinetic dynamics in the phase space. 

We shall here present only paradoxes that can be investigated from simple 
and mathematical analyzes, in view of that a large number of philosophical 
issues, less mathematical but not less important, have been debated rather 
intensively. 



2 A brief survey of the Boltzmann equation 

There are many descriptions of the Boltzmann equation in the literature [12- 
14]. For purposes of this paper, we shall recall the Boltzmann equation and its 
derivation in a brief and illustrative way. 

The core of the Boltzmann equation, which also serves as the basics of other 
kinetic equations, says that if collisions of molecules of a gas can be neglected in 
a certain region, the distribution function there satisfies the collisionless Boltz- 
mann equation 

df df F df 

at or m av 

which is sometimes called the continuity equation in the phase space. The proof 
of this continuity equation is obtained by considering a six-dimensional phase 
volume element (box-shaped) and formulating the total flux through all the 
facets of the volume element. According to textbooks, this equation is also 
equivalent to the path invariance of distribution function 



df 

dt 



= or /| r (t) lV (t) = Constant, (2) 

•m,v(t) 



where the subindexes r(t) and v(t) indicate that that the expression is valid 
along a molecule's trajectory. Conceptually speaking, equation (2) can be in- 
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terpreted as saying that such a gas is 'incompressible' in the phase space. The 
direct proof of the path invariance expressed by (2) is usually associated with 
a Jacobian formulation. The proofs of (1) and (2) in textbooks appear to be 
stringent except that there is one tacit assumption: the distribution function 
is perfectly continuous (differentiable, to be more exact) in terms of the six 
position and velocity components. Behind this continuity assumption, there is 
a more fundamental concept called molecular chaos [14]. It is explicitly or im- 
plicitly assumed that with help of molecular chaos the perfect continuity of all 
interested distribution functions is virtually ensured. 

To incorporate collisional effects into the formalism, collisions between two 
molecular beams, with velocities vi and v 2 respectively, are first considered. 
The number of type-1 molecules (with the velocity vi) leaving drdvi during a 
time interval dt is then found to be, according to the standard approach, 

dtdrdvi / u/(vi)/(v 2 )a-(vi,v 2 -> v' 1 ,v' 2 )dv2dv[dv 2 , (3) 

where and v 2 are the final velocities of type-1 and type-2 molecules respec- 
tively, u = |v 2 — v' x | = |v 2 — vi| and cr(vi,v 2 — > vi,v 2 ) is the cross section 
in the laboratory frame; likewise, the number of the molecules entering drdv\ 
during dt to be 

dtdrdv 1 u/(vi)/(v 2 )cr(vi, v 2 -> vi,v 2 )dv 2 rfvirfv 2 . (4) 

In view of the conservation of molecular number, the standard approach assumes 
that 

f a(Cl)dQ = f f <t(vi,v 2 -► v'^v'^dv'^v'z, (5) 

J £1 J J v' 2 

where cr(f2) is the cross section defined in the center-of-mass frame. The two 
cross sections in (3) and (4) are considered to be equal to each other because of 
the time-reversibility of collision. On this understanding, the net change of the 
molecular number per unit phase volume and per unit time can be expressed by 

/ «[/(vi)/(v 2 ) - /( Vl )/(v 2 )]a(fi)dnrfv 2 . (6) 

By regarding (6) as a correction term of the collisionless Boltzmann equation, 
the standard Boltzmann equation is obtained as 

S" +Vl '^ + ^- Jv7 = / Jf«)f( V 2) f(vi)f(v 2 )}ua(n)dQdv 2 . (7) 

The Boltzmann equation (7) is usually recognized as an independent and 
complete dynamical description of a gas consisting of hard-sphere molecules. 



3 



The terms 'independent' and 'complete' mean that the solution for the equa- 
tion can be constructed without getting help from other physical laws, such as 
Newton's equations for individual molecules. It is widely believed that if suffi- 
ciently powerful computers were available now, we would be able to solve the 
equation by applying certain types of finite difference schemes (under certain 
boundary conditions and initial conditions, of course). 

Although the Boltzmann equation and its derivations in textbooks seem 
quite stringent and have been accepted unanimously, many paradoxes related 
to it can be raised. Ironically enough, even a simple comparison between the 
left and right sides provides enlightening things to ponder. On the left side, 
molecules entering and leaving a six-dimensional phase volume element drdvx 
have been considered and there is a good symmetry between the position vector 
r and the velocity vector vi in terms of differential operations. On the right 
side, there is no symmetry between r and vi, all the integral operations are 
performed in the velocity space and the position vector r serves only as an 
inactive parameter. That is to say, although the derivation of the collisional 
operator is formally carried out in the phase space, only the velocity subspace 
has been in one's mind. If a forum were held in this community to discuss 
which space should be used in terms of formulating the collisional operator, a 
50-to-50 situation would be there: about a half of us may prefer the phase space 
in view of that the left side of the Boltzmann equation is indeed written down 
in the phase space, while the rest may choose the velocity space in view of that 
collision issues have, almost always, been investigated in the velocity space only. 
In Sect. 4 of this paper, we shall return to this subject. Much to our surprise, 
it will be unveiled that, unless some of our basic concepts are changed, we shall 
have intrinsic difficulties in formulating collisional effects. 

3 The continuity aspect of distribution function 

Here, we wish to show that the perfect continuity of distribution function is just 
an assumption; in many important situations, using certain types of discontin- 
uous distribution functions is absolutely a must. Throughout this section, in 
order to make the discussion less cumbersome, molecule-to-molecule collisions 
are considered negligible, at least in terms of the zeroth-order approximation. 

We first establish the connection between the perfect continuity of distri- 
bution function and the path invariance expressed by (2). Though the six- 
dimensional phase space is somewhat abstract, we can visualize the path in- 
variance in it quite easily. Fig. 1 is drawn to show the situation in the two- 
dimensional x — v x phase space; but things taking place in it can schematically 
manifest what happens in the full six-dimensional phase space. Suppose that 
the molecules inside the rectangle in Fig. la distribute rather uniformly at t = 0; 
these molecules will, due to molecular motion, distribute inside the parallelo- 
gram in Fig. lb rather uniformly at t = T. By excluding all external forces (to 
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make our understanding as simple as possible), we readily see that the area of 
the rectangle is equal to that of the parallelogram, which means that the molec- 
ular density inside the moving volume element keeps invariant. Based on the 
picture outlined above, we can use a sufficiently small phase volume element and 
let it move along a molecule's path, the average density of molecules in it will 
not change, irrespective of the shape and size of the employed volume element. 



Av x 



Ax 




(a) (b) 

Figure 1: Schematic of molecular motion in the phase space, (a) 
A number of molecules distribute inside the rectangle at t = 0. 
(b) They will distribute inside a paralellogram at t = T. 



Av x 



Ax 




(a) 



(b) 



Figure 2: For confined molecules, the situation can be much dif- 
ferent from that in Fig. 1. (a) All molecules distribute along a line 
at t = 0. (b) Separations between the molecules become larger at 
t = T. 

However, if the molecules are confined to a certain region in the phase space, 
the conclusion given in the last paragraph will be valid no longer. For this to 
be seen, investigate Fig. 2a, in which all the molecules distribute along one 
diagonal of the rectangle at t = 0, and they will, at t = T, move to new places 
and distribute along the diagonal of the parallelogram, as shown in Fig. 2b. If 
we have a phase volume element of a fixed size and shape and let it move with 
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the molecules, we find that the molecular density in it, defined as the molecular 
number divided by the phase volume, changes with time. Furthermore, it can 
be seen that if the shape and size of the employed volume element are allowed 
to change, we shall get different molecular densities, varying drastically from 
zero to infinity. 

A crucial question arises. Can molecules in a real gas be so confined? To find 
out the answer, consider a beam of molecules striking a surface S and 'reflected' 
from it; or consider a surface S moving toward a group of stationary molecules 
and driving them to move. As shown in Fig. 3a, if S is convex and relatively 
smooth, the molecules will, after colliding with the surface, spread out in space 
with definite velocities. (Even if the surface is not truly smooth, the picture 
shown will not change significantly.) By counting the molecules in a moving 
spatial volume element Ar and in a velocity range between v and v + Av, 
we find that the molecular density within ArAv decreases. This means that 
these molecules are diverging in the phase space or, in terms of the distribution 
function / describing them, we get an expression as 



dt 



r(t),v(t) 



df df F df 
<0 or -± + v. 1 ± + -.-L <0. 
ot or m ov 



(8) 



For the situation shown in Fig. 3b, we may have, in a certain region along 
molecular paths, 

df >0. (9) 

r(t),v(t) 



dt 



Equations (8)- (9) and Fig. 3 tell us that diverging and converging molecules asso- 
ciated with discontinuous distributions can be constantly produced by collisions 
between molecules and solid boundaries and these discontinuous distributions 
will keep discontinuous along the molecular paths. 



(AS); 




S/ 



(AS)} 



(a) 



(b) 



Figure 3: Illustration of how boundaries 'reflect' incident 
molecules, (a) If the surface is convex and rather smooth, the 
reflected molecules diverge, (b) The reflected molecules may con- 
verge. 
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Are the diverging and converging molecules discussed in the last paragraph 
important to gas dynamics? Of course, they are. If the gas in question is a really 
rarefied one, almost all the molecules in it can be regarded as the ones produced 
directly by certain boundaries, which means that before a good equilibrium is 
established a great number of molecules need to be described by discontinuous 
distributions. For a more practical gas, we can at least say that the regions 
immediately near boundaries are very much filled with diverging and converging 
molecules abnormal to the Boltzmann equation, which is sufficient to challenge 
the standard methodologies in view of that the distributions in the inner region 
(assumed to exactly obey the Boltzmann equation for the discussion herein) 
has to be connected with the distributions caused by boundaries in one way or 
another. 

We continue searching for cases in that distribution functions are not per- 
fectly continuous. In Fig. 4a, on the left side of the wall the space is filled with 
a dense gas and on the right side with a dilute gas (or just a vacuum). It is 
obvious that the molecules leaking through a small hole on the wall will behave 
like diverging molecules in the phase space. In Fig. 4b, a gas with a very high 
density is confined to a small spherical ball. Suppose that the ball confinement 
suddenly vanishes at a moment (chosen as the initial time), and then the gas 
expands almost freely. Obviously, these expanding molecules also behave like 
diverging molecules. Finally, referring to Fig. 4c, consider two narrow beams 
of molecules meeting with each other in a small region. It can be seen that 
the scattered molecules by the collisions there will diverge in the phase space. 
The examples presented thus far tell us that molecule-to-boundary collisions, 
molcculc-to-molcculc collisions, boundary conditions and initial conditions are 
all capable of making molecules diverge or converge in the phase space. 




(a) (b) (c) 

Figure 4: (a) A gas leaking out of a small hole, (b) A gas expands 
freely after a certain time, (c) Two molecular beams collide in a 
small region. 

As far as the discontinuity issue is concerned, it is quite instructive to inves- 
tigate the path invariance theorem (2) more thoroughly. Though the theorem 
is intimately related to the continuity assumption of distribution function, it 
can, in a subtle and interesting way, lead us to betray the existing theory and 
admit the necessity of using discontinuous distribution functions. To find out 



7 



true implication of the path invariance, we rewrite it as, with all external forces 
neglected for simplicity: 

/(t,r,v) =/[to,r ,v )] =/[to,r-v(t-t ),v], (10) 

where to is regarded as the initial time (to < t). Based on equation (10), we 
may call r = r — v(t — t ) a source point and r an image point. We then assume 
the distribution function at to to be known and try to determine the velocity 
distribution at a specific point r and at a later time t, denoted by ft, T (v,fl) 
where v = |v| and the solid angle specifies the direction of v. From equation 
(10), it is seen that for a certain value of v all the source points of f t , T (v,Fl) 
distribute on the spherical shell that has the radius v(t — to) and is centered at r. 
Since v is allowed to vary from zero to infinity, all points in the spatial space can 
serve as the source points of /t jr (i>, fi). In other words, the path invariance says 
that the velocity distribution at any spatial point is a condensed image of the 
entire spatial space; or, to say it in a converse way, any local spatial structure 
will constantly send its image to participate in forming the velocity distribution 
at every and each point elsewhere. A direct, and somewhat surprising, corollary 
of the above statements is that if a local structure is discontinuous no point in 
the space is truly immune from discontinuity. 




(a) (b) 

Figure 5: At each point in space, the velocity distribution is 
formed as an image of structures elsewhere. The stuctures can 
be (a) a solid block surrounded by a gas; (b) a continuous gas 
with density gradient. 

Making things worse, even a continuous distribution can constantly send 
'bad images' to other places. Let's compare between Fig. 5a and Fig. 5b. In 
Fig. 5a, there is a solid block surrounded by a continuous gas. The velocity 
distributions at other points, as images of the block and the gas, must involve 
corresponding discontinuity, as analyzed in the last paragraph. We then turn 
to Fig. 5b and assume that the distribution function in the oval region can be 
expressed by a local Maxwellian n(r) cxp(—mv 2 /2T), where the temperature T 
is constant and the symbol n(r) represents that the density there has a certain 
spatial gradient (but still perfectly contitinuous). From (10), we find that the 
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sharpness of the velocity distribution at an image point r outside the region can 
be characterized by 



df(y) 



n(r 2 ) - n (ri) 
AO 



-e 



2 /(2T) 



(11) 



where ri and r 2 are two source points in the oval region and is the angle formed 
by the two paths r — ri and r — r 2 (both are straight lines with no external force). 
This expression informs us that as the image point r goes farther and farther 
away from the region (while holding ri and r 2 fixed), the profile of velocity 
distribution will have a sharper and sharper gradient in terms of the velocity 
angles and, finally, get almost the same shape as that caused discontinuous 
source structures. For convenience of discussion, we shall say that the veloc- 
ity distributions illustrated in Fig. 5b are those involving quasi-discontinuity, 
though similar phenomena have been studied in mathematics under the title 
'nonuniform continuity' [15]. 

The last two paragraphs manifest several interesting characteristics of gas 
dynamics, (i) Macroscopic structures constantly send their images (structure 
information) to shape and form microscopic structures elsewhere, (ii) When 
the structure information is conveyed through space it is condensed and sharp- 
ened progressively, (iii) Since a local distribution is a condensed image of its 
surrounding environment, it is improper and unnecessary to endow distribution 
functions with special continuity properties, (iv) Physical events make their in- 
fluences through one-dimensional trajectories, which is described by collisionless 
Newton's equations or Hamilton's equations. 

We now wish to examine the discontinuous distribution function produced 
by a boundary of finite size in an analytical manner. Experimental facts tell 
us that collisions between molecules and boundaries cannot be truly elastic; 
stochastic and dissipative forces must get involved to a certain degree[l][16]. To 
describe molecules 'emitted' by such boundaries, it is adequate, in particular 
for a rarefied gas, to define the local emission rate a in such a way that 

adtdSdvdfl (12) 

represents the number of molecules that are reflected by the surface element dS 
and emerge in the speed range dv and in the solid angle range dQ during the 
time interval dt. Referring Fig. 3 and dividing S there into N elements, we find 
that for the ith element of S the reflected molecules are like ones emitted from a 
point source and the corresponding distribution function at a point r elsewhere 
can be expressed by, with external forces neglected, 

/ (r, v, fi) = a{ASr j; 3 Vi(v)8(n - O r _ rot ), (13) 
|r - roiftr 

where f2 stands for a solid angle in the velocity space at r, Vi(v) is a function of v, 
r i is the position vector of (AS), and the solid angle fl r -r oi takes the direction 



9 



in that a reflected molecule will continue to move from r after passing through its 
trajectory roi — > r. Since external forces are ignored, f2 r _ r(H simply points in the 
direction of r — tqi. In terms of r and v, this distribution function is continuous; 
in terms of fl, it is confined to a certain 'point' f2 r _ r(H . In this sense, we wish 
to say that the distribution function (13) is a function of four dimensions. The 
distribution function at r associated with all reflected molecules from all the 
area elements is 

/(r, v, Q) = , a{AS ^ 3 Vi(v)6(Q O r _ ro J. (14) 

Note that, on the usual understanding of distribution function, the summation 
symbol in (14) is a formal one and we cannot replace it with the integral sign 
since each of ^ r _ roi represents a distinctive direction. Fig. 3 has shown that 
the molecules related to (14) can diverge or converge in the phase space, mainly 
depending on how S is arranged. As a limiting case, if S is sufficiently small, 
equation (14), being identical to (13), represents diverging molecules. 

Though the above discussion appears to be trivial, many new concepts can 
arise from it. As one thing, allowing ^-functions to express distribution functions 
shakes the prevailing belief that distribution functions must be continuous and 
gas dynamics must be analyzed by differential equations. As another, with the 
use of such (5-functions, path information of molecules (which is described by 
collisionless Newton's equations or Hamilton's equations) enters into the picture 
in a natural and treatable way. Finally, in order to relate this type of distribution 
function to the usual type of distribution function, we are actually compelled 
to invoke integral-type formulations. 

Before finishing this subject, one issue may come into our mind. In fluid 
mechanics, several types of discontinuity, such as those related to shock waves, 
have been treated by finite difference schemes in a fairly successful way, why 
should we give so much concern to discontinuity in the kinetic regime? There are 
a number of things that can justify our attitude. Firstly, unlike what happens 
in fluid mechanics, discontinuities discussed in this paper spread out and involve 
almost all points in space. Secondly, we have seen that the quasi-discontinuity 
produced by distribution functions will become sharper and sharper when prop- 
agating, and this type of dynamics cannot be dealt with by finite difference 
schemes. Finally, all types of discontinuity presented in this section transfer 
along molecular paths, suggesting that path information must play a vital role 
in the would-be formalism, whereas path information is disregarded in almost all 
finite difference schemes of solving the Boltzmann equation (see the discussion 
in Sect. 2 also). 
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4 The collision aspect of the Boltzmann equa- 
tion 



The collisional operator of the Boltzmann equation seems to have been 'clearly 
elaborated' in many textbooks. However, as will be seen, these elaborations 
actually involve vague descriptions and tacit assumptions. If we look at them 
with a mind less preoccupied, paradoxical concepts hidden in them will unveil 
themselves. 




Figure 6: Schematic of the time reversibility of two individual 
molecules, (a) The original collision; and (b) the inverse collision. 



To begin with, we recall the usual description of time reversal for a collision 
of two molecules with the same mass (still distinguishable according to classical 
mechanics). The initial and final velocities of the first molecule are denoted by 
vi and v[ respectively; the initial and final velocities of the second molecule by 
v 2 and v 2 respectively. The following two collisions, referring to Fig. 6, 

Vi,v 2 ->vi,V2 and - v[, -v 2 -> -Vi, -v 2 (15) 

can be deemed as a pair of corresponding collisions. Under this understanding, 
the time reversibility of collision implies that if the first collision of the pair is 
physically possible then the second one must also be physically possible. This 
type of time reversibility has been solidly established in classical mechanics and 
we shall discuss it no more. 

Now, we try to investigate whether or not there is time reversibility in terms 
of collisions between two molecular beams. Apparently, to the Boltzmann equa- 
tion a time reversibility of this type is much more relevant. In what follows, 
it will be shown that, much contrary to the conventional wisdom, no time re- 
versibility of this type can be found out: neither an intuitive one, nor a mathe- 
matical one. 

Intuitively, we may consider the following pair of corresponding pictures as 
a 'candidate' for such time reversibility. The original collision picture is that 
two molecular beams with definite velocities collide with each other and the 
molecules scattered by collisions diverge in the spatial space with many differ- 
ent velocities, as shown in Fig. 7a. The inverse collision picture is that many 
different converging beams collide with each other and the molecules scattered 
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by collisions form two molecular beams with two definite velocities, as shown in 
Fig. 7b. In no need of much discussion, we all know that the first picture makes 
sense in statistical mechanics, while the second one does not. 




(a) (b) 

Figure 7: A possible candidate of the time reversibility of beam-to- 
beam collision: (a) the original collisions; and (b) possible inverse 
collisions. 

Mathematically, according to the textbook treatment [12], the time reversibil- 
ity of beam-to-beam collision can be expressed by 

cr(vi, v 2 -» vi, Va) = cr(vi,v 2 -> vi,v 2 ), (16) 

where the cross section <r(vi,v 2 — > v' 1; v 2 ) is defined in such a way that, after 
collisions between a beam of type-1 molecules at vi and a type-2 molecule at 
v 2 , 

N = cr(vi,v 2 -» v^v^rfvirfvi (17) 

represents the number of type-1 molecules emerging between and v[ + dv[ 
per unit incident flux and unit time, while the type-2 molecule emerges between 
v 2 and v 2 +dv 2 ; and the cross section cr(v[, v 2 — > Vi, v 2 ) is defined in the same 
way. 

An unfortunate fact is that the concepts related to (16) and (17) suffer from 
fatal problems. For the collisions, the energy and momentum conservation laws 
state that 

v i + v 2 = v i + v 2 = 2c and |v 2 - v^| = |v 2 - Vi| = u, (18) 

where c is the velocity of the center-of-mass and u is the relative speed of 
the two molecules. Fig. 8a shows how we can get c and u from the initial 
velocities vi and v 2 , while Fig. 8b shows how c and u form constraints on the 
final velocities v' x and v 2 . It is clear that the final velocities of all scattered 
molecules must fall on a spherical shell S of diameter u = |u| in the velocity 
space, which will be called the accessible shell. Keeping this accessible shell in 
mind, two misconcepts connected with the definition (17) manifest themselves. 
The first is that after dv[ is specified, specifying dw' 2 in the definition (17) is a 
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work overdone. The second is that the cross section should not be defined in 
reference to a velocity volume element (like dv[ in the definition) . If we insist on 
doing so, the resultant cross section will unpleasantly vary from zero to infinity, 
depending on the shape and size of the employed velocity volume element. 




(a) (b) 

Figure 8: Constraints imposed by the energy and momentum con- 
servation laws, (a) The velocity c and u — |u| are determined 
by the initial velocities, (b) The final velocities of the scattered 
molecules and v' 2 must fall on the accessible shell S of diameter 
u. 



■e-, s 



dv[ 



dv[ 



S 



(a) (b) (c) 

Figure 9: Possible velocity volume elements dv[ are given, each 
of which encloses a small part of the accessible surface S. They 
can be (a) a spherical ball, (b) a slim cylinder, and (c) a tabular 
cylinder. 

To see the conclusion given in the last paragraph more vividly, look at Fig. 9, 
in which several possible shapes have been given to dv[ . For the situation shown 
in Fig. 9a, the number of molecules entering the spherical ball, chosen to be 
dv[, per unit time can be expressed by N « pnr 2 , where p is the area density 
of molecules on the accessible shell caused by unit flux of type-1 molecules and 
r is the radius of the spherical ball dv[. By noting that p must be finite, we 
find that the cross section defined by (17) is equal to, with dv' 2 neglected, 

_iV_ _ pirr 2 3p 
a dv[ 47rr3/3 4r' { 1 
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and it tends to infinity as r — > 0. Nevertheless, if the volume element dv[ is 
chosen to be the one in Fig. 9b, then a tends to zero; if the volume element dv[ 
is chosen to be the one in Fig. 9c, a tends to infinity again. These examples 
simply mean that the cross section defined by (17) is, in fact, an ill-defined 
quantity, and thus the time reversibility (16) defined by such cross sections is 
also groundless. 

It is now clear that the Boltzmann collisional operator indeed involves mis- 
concepts. An interesting and essential question arises: Can we still formulate 
collisional effects in the context of the standard theory? In the rest of this 
section we try to get an answer to this question by reinvestigating collisional 
effects in the three-dimensional velocity space and in the six-dimensional phase 
space. Surprisingly enough, such investigations will give us nothing but more 
paradoxes. 

Let's first investigate the number of the molecules that enter a volume ele- 
ment in a velocity space, denoted by dvi herein. This type of investigation is 
not truly relevant, since the Boltzmann equation is supposed to be cast in the 
full six-dimensional phase space, but we take this approach anyway, in order 
to see where we arrive at. Following the usual methodology, we consider two 
molecular beams with velocities and v 2 (called type-1 and type-2 molecules 
respectively), and examine how many type-1 molecules emerge, after collisions, 
in the range between Vi and vi + dv\. As the discussion after (18) has already 
manifested, because of the energy and momentum conservation laws, all the 
type-1 molecules scattered by collisions must fall on the accessible shell. If we 
insist on knowing the density of the molecules entering the volume element dv\ , 
we get a result that varies drastically from zero to infinity, depending on the 
position, size and shape of dv\. 

We now turn our attention from the velocity space to the phase space and try 
to examine the net change of the molecular number in a six-dimensional phase 
volume element due to collisions. For the reason that has been clarified, the 
following examination will be done without getting help from the cross section 
<t(vi,v 2 -> vi,v 2 ). 

Firstly, the molecules leaving drdv\ during dt due to collisions are of our 
concern. Rather conventionally, we identify /(vi)dvi and /(v2)c?V2 as two 
molecular beams. By switching to the center-of-mass frame and studying the 
collisions related to the two beams, the number of the collisions in dr during dt 
can be directly obtained as 



where <r(Q) is well-defined in the center-of-mass frame[17]. Integrating the above 
formula over v 2 yields the total number of the collisions taking place in drdvi 
during dt 




(20) 




(21) 
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Until now, every step is almost the same as that in the standard approach, and 
the comparison between this expression and the Boltzmann collisional operator 
seems to state that things go and will go as smoothly as expected. However, a 
careful inspection will tell us different stories. The form of dtdrdvi in expres- 
sion (21) reminds us that dt, dr and dv\ are three differentials; in obtaining 
the Boltzmann equation, each of them must be allowed to approach zero inde- 
pendently. Conflicting conclusions will be drawn when different ways of taking 
these limits are adopted. One of the conclusions is that: if \dr\ » \v\dt\, 
where |c?r| stands for the length scale of dr, almost all the molecules not suffer- 
ing collisions will stay inside drdvi during dt and the molecules expressed by 
(21) are indeed the ones that get velocity changes due to collisions and leave 
drdv\ during dt. But, another conclusion is quite opposite: if |c?r| << |vidt|, all 
the molecules initially inside drdvi will get out of drdvi in a time much shorter 
than dt, irrespective of suffering collisions or not (see Fig. 10a). In the latter 
case, it is meaningless to say that the molecules described by (21) are those 
leaving the volume element during dt just due to collisions. In other words, to 
make (21) relevant, we have no choice but to assume that |dr| >> |vi<fi|. An 
unfortunate fact is that no sound reason can be found out for that we can prefer 
this assumption to its converse. 




Figure 10: Schematic of what happens in an infinitesimal element 
drdvi during dt. (a) If |dr| << |vidt| all molecules will leave dr 
almost instantly, (b) If |dr| << |vidt| most of molecules produced 
by collisions will not stay within dr after dt. 

Then, we examine molecules entering a phase volume element drdvi during 
dt due to collisions. Following the textbook treatment, two beams with velocities 
v[ and v' 2 are assumed to collide within a spatial volume element dr during dt, 
giving contribution to the molecules expressed by f(r, v\,t)drdv\. Since no time 
reversibility can be utilized as a 'shortcut', we have to try the task from scratch. 
It turns out that not one but many paradoxes will be in front of us. As the first 
paradox, if we assume, by adopting the customary thought, that only collisions 
taking place within dr are worth considering, then the size and shape of dr will 
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become irrelevant and the six-dimensional collision issue will be reduced to the 
three-dimensional collision issue in the velocity space. As have already been 
discussed, the number of the molecules entering a three-dimensional velocity 
volume element is an ill-defined quantity and we have no way to determine it 
definitely. As the second paradox, even if we were able to formulate the relevant 
molecular number produced by collisions within drdvi during dt, we would still 
encounter serious difficulties. For instance, we again need to let dt, dr and 
rfvi approach zero. If |dr| << |vidt|, all produced molecules will 'instantly' 
leave dvdvi (see Fig. 10b) and cannot be treated as a true contribution to 
/(r,vi,t)drdvi. As the third paradox, molecules produced by local collisions 
must diverge in the phase space, as revealed in the last section; such molecules 
have to be treated with care and cannot be regarded as a regular contribution 
to /(r, vi,i)drdvi. 

The discussions in this section, as well as some of those in the last section, 
suggest one essential and delicate thing: the phase space is a space much more 
complex than it appears and the collective behavior of molecules in it has to be 
studied in a careful and comprehensive way. The phase space, in which x,y,z 
and p x ,Py,Pz are regarded as variables independent of each other, has been em- 
ployed for many years and has served as a good mental and mathematical tool in 
classical mechanics. Its success is actually related to the existence of Hamilton's 
equations. With help of Hamilton's equations, the true correlation between t, r 
and v will be restored when the solution of Hamilton's equations is obtained. In 
statistical mechanics, the situation is subtly different. As has been illustrated 
in this section, to formulate collisional dynamics, it is sometimes necessary to 
consider the correlation between t, r and v and sometimes not necessary to do 
so. Before and after collisions, the molecules obey Hamilton's equations for 
individual molecules, in which dt, dr and dv indeed correlate. When collisions 
take place, the molecules are governed by laws of beam-to-beam collisions, in 
which certain types of statistical assumptions have to be adopted and Hamil- 
ton's equations for individual molecules do not play an essential role. Though 
the standard collisional operator tries hard to combine these two mechanisms, 
some conflicts remain there. 

5 Conclusions 

In summary, two paradoxical aspects of the standard Boltzmann equation have 
been presented. The first paradoxical aspect is related to the usual concept 
of distribution function. An tacit assumption of the existing kinetic equations 
is that distribution functions, though describing discrete molecules, must be 
mentally and practically continuous. The discussions in this paper, however, 
show that distribution functions of realistic gases can have complex structures 
over infinitesimal range, and such distribution functions do not keep invariant 
along a molecule's path even if the local collisions can be ignored. The second 
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paradoxical aspect is related to the standard collisional operator. It is pointed 
out that the time reversibility in terms of bcam-to-beam collisions is actually 
fictitious. Along this line, it is unveiled that there are unsurmountable diffi- 
culties in terms of formulating molecules entering and leaving an infinitesimal 
phase volume element. All these things pinpoint one thing in common: it is not 
possible to formulate collective behavior of molecules in the phase space mainly 
based on a differential-type analysis. 

Many more fundamental questions arise sharply. Some of our recent works 
try to make more analyses and put forward alternative approaches[18, 19]. With 
help of a development in quantum mechanics [20], part of the discussion in this 
paper can be extended to the regime of quantum statistical physics. 

This paper is supported by School of Science, BUAA, PRC and by Education 
Ministry, PRC. 
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